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Abstract 



We have studied the discrete nonlinear Schrodinger equation (DNLSE) 
with on-site defects and periodic boundary conditions. When the array dy- 
namics becomes chaotic, the otherwise quasiperiodic amplitude correlations 
between the oscillators are destroyed via a symmetry breaking instability. 
However, we show that synchronization of symbolic information of suitable 
amplitude signals is possible in this hamiltonian system. 
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1. Introduction. The subject of our study is a subtle form of synchronization that occurs in 
the discrete nonlinear Schrodinger equation (DNLSE). The DNLSE describes, among other 
systems, localized modes in long protein systems , arrays of nonlinear mechanical pendula 
and the propagation of light in nonlinear optical waveguide arrays based on Kerr media 
[1]. These optical arrays were studied intensively in the past several years, mostly in the 
context of intrinsec localized modes also known as breathers, both theoretically [2,3] and 
experimentally [3]. The DNLSE also arises in recent models of Bose-Einstein condensates 
trapped in periodic potentials generated by optical standing- waves [4]. 

In hamiltonian dynamical systems, the usual notion of synchronization cannot be applied 
, since volume must be preserved in phase space [5,6]. Indeed, the trajectories of coupled 
systems have to collapse to the synchronization manifold [5,6], which is only possible in 
dissipative systems. However, a particular form of synchronization in hamiltonian systems 
was found in article [7], where the orbits of two coupled systems cover almost the same 
region of the individual phase space but , typically, at different times [7]. In the present 
work, we find that another form of synchronization is possible in a hamiltonian system such 
as the DNLSE. We refer to synchronization of symbolic information (SSI). This concept 
has been applied recently to describe the simultaneous production of common information 
by two coupled chaotic oscillators as diverse as a logistic map and an electronic circuit [8]. 
Our article has four parts. The DNLSE and a new set of DNLSE solutions are discussed 
in the second part. In part three, we consider the DNLSE with an on-site defect, where 
chaotic synchronization of symbolic information occurs. Finally, in part four, we give the 
conclusions. 

2. The Model and a Family of Stationary Solutions. We will consider the DNLSE within 
the framework of optical waveguides. Let us consider an array of one- dimensional coupled 
waveguides. The neighboring waveguides are separated from each other by the same distance 
d and , therefore, the coupling constant between these is the same. We consider the case 
with lossless Kerr media. Within the framework of the coupled mode theory, the evolution 
of E n , the electric field envelope in the nth waveguide, is given by the following equation 

2 



[2]: i^f + !3 n E n + C(E n -i + E n+1 ) + 7 | E% \ E n = 0. This system is assumed to have 
periodic boundary conditions. In this equation, f3 n is the linear propagation constant of 
the nth waveguide, C is the linear coupling coefficient and 7 is the nonlinear parameter. 
By introducing the dimensionless field Q n = ^J r y/2CE n exp(— i((5 + 2C)z), this equation 
transforms into the discrete nonlinear Schrodinger equation (DNLSE) given by : 

■ dQr, 



+ KQn + (Qn-l + Qn+1 ~ 2Q n ) + 2 | Q> n \ Q n = 0, (I) 



where ( = Cz is the normalized propagation distance and 5 n = ^ n ~^ , The hamiltonian 
equations related to the DNLSE can be obtained by transforming into action-angle variables 
(I n , On), where Q n = ^/T^exp(—i9 n ). The quantity I n stands for the light intensity of the 
nth waveguide [1]. 

Now, we will consider a family of stationary solutions for which 5 n = 0. To find these, 
we deal with the nonlinear map approach [9,10]. We underline that there are different ways 
to find stationary solutions of the DNLSE [1]. In article [9,10], it was shown that homoclinic 
and heteroclinic orbits of a suitable hamiltonian map support breather solutions in the 
DNLSE. Instead, we will consider those stationary solutions which are determined by the 
(elliptic) stable periodic orbits of the hamiltonian map mentioned above. These two types of 
DNLSE solutions are clearly different. Indeed, these elliptic and homoclinic orbits , as it is 
well known, have different origins and properties [11]. While the first are periodic, according 
to the Poincare-Birkhoff theorem, the homoclinic orbits have no periodicity, i.e., repeated 
iteration of the associated map produces new homoclinic points [11]. These new stationary 
solutions are obtained by making 4£k = Q and 9 n = 9 m for any n ^ m. Moreover, we can 
define the frequencies on these tori I n = P% by setting ^ = A, where A is a parameter. 
Hence, the following map is obtained : X n +i = P n ', P-n+i — (r n — 2 | P n \ 2 )P n — X n , 
where T n = 2 — A — 5 n . The Jacobian J of this map is area preserving , i.e. J = 1. 
For the sake of definition, we will refer to this map as the cubic map (CM). To find our 
stationary solutions, we will consider the case T = 2 — A for which 8 n — 0. This map has 
the symmetry X n — * — X n and P n — > — P n . The fixed points of the CM are : (0, 0) and 



(±y/T/2 - l,±y/T/2 - 1). When T = 2.5, we find a saddle point at (0,0) and two elliptic 



points at (±wT/2 — l,±wT/2 — 1). In Fig. la, quasiperiodic orbits surround the elliptic 



point (•y/r/2 - 1, ^r/2 - 1) and further away, we identify a period seven island chain and its 
resonances. Surrounding these resonances, we also observe the characteristic chaotic sea in 
this nonintegrable map [11]. We choose the stationary field P n associated with the resonance 
of period seven. This is shown in Fig. lb. The stability of P n was numerically studied by 
integrating the DNLSE. We have carried out these integrations with slightly different initial 
conditions at £ = 0. That is, | Q n (( = 0) |= P n + v n , where v n stands for a small random 
perturbation between — u max < u n < u max and | 9 n (0) |<C 1 for any n. In our calculations 
Vmax = 10~ 3 . We find that the standard deviation of is proportional to the initial 

deviation from the stationary solution. This is a result of the stability of this stationary 
solution. The finite number of power spectrum components of indicates that these 

oscillations are not chaotic. Indeed, we numerically found that all the Liapunov exponents 
A n of this system converge towards zero. The power spectrum S(F) of /3(C) is shown in the 
dashed line of Fig. 2a, where F is the renormalized frequency. 

3. Synchronization of Symbolic Information (SSI). To find the aforementioned correlation 
effects, we add a small defect to the linear propagation constant of the central waveguide, 
which in Fig. lb has precisely the largest amplitude, i.e. £3 7^ and 5 n = for n 7^ 3. 
We make use of the stationary solution P n of period seven as the initial condition for this 
perturbed DNLSE. To estimate the extent of coherence between the intensities of a given pair 
of waveguides (I m ,I n ) we consider the linear cross correlation function R of this waveguide 
pair. Typically, for a small enough positive defect < 5 3 C 1, the intensities are not chaotic 
and the symmetric waveguide pairs , i.e. (Ji,J 5 ), (hjh) and (Iq,Iy), become strongly 
correlated for a suitable 5 3 > 0, i.e. R ~ 1 in the three waveguide pairs [12]. 

However, in an interval of the negative defect £3 < , it is possible to induce hamiltonian 
chaos in this system. In contrast to the previous case, the onset of chaos is manifested when 
the intensity pairs (Ji, J 5 ), (hjh) an d {Iq,It) are no longer correlated. The exponential 
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divergence of nearby trajectories explains why two otherwise intensity correlated waveguides 
, such as ii and J 5 , lose their relative symmetry when chaos arises. This divergence is shown 
in Fig. 2b, in which each of the aforementioned intensity pairs executes initially stable small 
amplitude oscillations before displaying erratic behavior. This behaviour is similar to that 
near the coupling resonance in the three-dimensional billiards problem, where a particle 
bounces back and forth between a smooth and a periodically rippled wall [11]. 

In the chaotic regime, in spite of the lack of symmetry of the aforementioned intensity 
pairs, a very interesting form of synchronization takes place between certain signals emerging 
from different pairs of waveguides. To this end, let us define the variables : J15 — Ii — I 5 , 
J24 = h — h an d J 76 = I-j — Iq. Our simulations show that J15, J 2 4 and J^ 6 have the same sign 
during most of the time. In fact, the signs of J mn are different during a negligible interval. 
This is clearly appreciated in Fig. 2c. We would like to know if the collective behaviour 
of the signals J mn suggests some kind of synchronization. The answer is affirmative from 
the point of view of synchronization of symbolic information [8]. According to this notion, 
two arbitrary oscillators are perfectly synchronized ,in an information sense, if they produce 
the same information, i.e., symbols generated by one system map one-to-one to symbols 
emitted by the other system. Strictly speaking, this form of synchronization requires that 
the common information be emitted at precisely the same time. In article [8], the chaotic 
signals of two coupled systems are compared using their symbolic dynamics. In our case, 
Fig. 2c suggests that the signals Ji 5 , J 2 4 and J 76 exhibit equivalent information at the same 
rate. This ,however, does not contradict the fact that the usual notion of synchronization 
[5,6] cannot be applied to hamiltonian systems ,such as ours. Indeed, in our system, we 
compare the chaotic signals J nm using just their symbolic dynamics. This is along the lines 
of the notion of event synchronization [13], where the relative timings of events in the time 
series of different physiological signals are compared. 

The signals Ij, where j = 1,..,7, are chaotic as suggested by the broadband power 
spectrum S(F) of /3(C) shown in the solid line of Fig. 2a. To generate a symbolic sequence 
out from Ji 5 , J 2 4 and Jjq, we replace the value of by " — 1" or "1" if Jy < or > 
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respectively. Let us call these signals Kij. In addition, the symmetry of the system indicates 
that < Jij >=< >= 0, where <> stands for the sample average. We have calculated 
the autocorrelation function (ACF) C(S) for each of these signals, where S is the space lag 
[11]. The length of these signals is roughly C ~ 10 5 . The ACF C{S) of K 15 , K 24 and K 76 
coincide exceptionally well for all space lags S considered. Indeed, all of them collapse to 
the solid line of Fig. 2d, where £3 = —0.0078. In this figure, the same result was obtained 
for 5s = —0.0158. In both cases, the associated linear cross correlation function R for any 
pair of signals K mn is given by R 1. In contrast, the ACF C(S) of J 15 , J 2 a and J 76 agree 
only at a qualitative level. Moreover, the first zero of these ACF C(S) occurs at the same 
lag S, since these series change sign almost simultaneously as shown in Fig. 2c. 

It is remarkable how the chaotic synchronization of symbolic information (SSI) of J 15 , 
J 24 and J 76 is manifested in the dynamics of A9 mn = m — 9 n when 5 3 < 0. Indeed, as shown 
in Fig. 3a, all the A9 mn are bounded when SSI occurs , i.e. when R « 1 for any pair of series 
K mn . This locking of the phases 9 m and 9 n , ] 9 m — 9 n \< 2tt when | 9 m (0) — 9 n (0) |< 2tt, 
is also referred to as phase synchronization [14]. The latter has been extensively studied 
in the context of coupled self-sustained chaotic oscillators [6,14]. Moreover, when 5 3 < 
and the onset of phase slips of A9 mn occurs, degradation of SSI (0 < R < 1) takes place. 
This is observed in Fig. 3b and Fig. 3c for which £3 = —0.25 and £3 = —1.5 respectively. 
When the dynamics of A9 mn looks like a random walk, SSI is no longer a typical feature 
of the series J mn . This is clearly appreciated in Fig. 4a, where the trajectory , at discrete 
intervals of (, is projected on the plane (Jis,^)- If SSI took place, the dots of this plane 
were only in the first and third quadrant of the coordinate plane as shown by Fig. 4b and 
Fig. 4c. To further characterize SSI , let us consider the spectrum of Liapunov exponents 
A n . This is shown in Fig. 4d. The DNLSE has two constants of motion, namely the norm 
and the hamiltonian and , moreover, this flow is autonomous. This, along with the symmetry 
of the Liapunov exponents of hamiltonian systems [11], suggests that at least four Liapunov 
exponents are equal to zero. For 5 3 = —0.0078, there is a single positive Liapunov exponent 
whose magnitude is much larger than that of the other positive exponents. However, for 
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#3 = —0.25 and £3 = —1.5 five positive Liapunov exponents have roughly the same order of 
magnitude as shown in Fig. 4d. Finally, in the context of the aforementioned BEC models, 
the phase locking of A9 mn is associated with the superfluid regime. However, when A9 mn is 
unbounded, the system is said to behave as an insulator [4] 

4- Conclusion. We have studied a subtle form of chaotic synchronization that arises in a 
hamiltonian dynamical system with several degrees of freedom. This behaviour emerges in 
a family of solutions of the DNLSE, which have as initial conditions the neighborhood of 
certain stationary solutions. The latter, in turn, are given by the resonances of a suitable 
hamiltonian map. In a small negative interval of the defect (63 < and | £3 |<C 1), suitable 
chaotic signals, which are generated by different pairs of oscillators, synchronize in the infor- 
mation sense. Indeed, almost the same binary symbolic dynamics is emitted simultaneously 
by these signals. When this happens, phase synchronization of all oscillators follows and 
, roughly, just a single positive Liapunov exponent is present. For arbitrary initial conditions 
or parameters, typically, synchronization of symbolic information does not take place. It is 
possible to extend this approach to lattices with periodic defects and cuadratic nonlinear- 
ities. In view of the wide applicability of the DNLSE, we believe that the aforementioned 
synchronization effect is generic in lattice models, 
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FIGURES 

FIG. 1. (a) Plot of the stationary field P n versus P n +i for T = 2.5. (b) Plot of the stationary 
field amplitudes P n versus waveguide index n for the period seven resonance and T = 2.5. 

FIG. 2. (a) Plot of log 10 S(F) versus log 10 (F) for I 3 ((). v max = 10" 3 , 6 3 = -0.0078 (solid 
line) and 5 3 = (dashed line), (b) Plot of the intensities I\ (solid line) and 1$ (dashed line) versus 
£, where £3 = —0.0078. (c) Plot of J15 (solid line) J24 (dashed line) J-j% (dotted line) versus C for 
5 3 = -0.0078. (d) Plot of the ACF C(S) of K 15 (() for <5 3 = -0.0078 (solid line) and <5 3 = -0.0158 
(dashed line). The function C(S) of Kx§, K24 and K-jq coincide exceptionally well. 

FIG. 3. (a) Plot of A0i5, A6 2 4, and A6 76 versus ( for <5 3 = -0.0078. Plot of A0 15 (solid line), 
A024 (dashed line) , and A6jq (dotted line) versus ( for (b) S 3 = —0.25 and (c) 6 3 = —1.5. 

FIG. 4. Plot of J15 versus J 24 and J 76 for (a) <5 3 = -1.5, (b) 5 3 = -0.0078 and (c) 5 3 = -0.0158. 
(d) Spectrum of Liapunov exponents A n for <5 3 = —0.0078 (solid line), S 3 = —0.25 (dashed line), 
and e>3 = —1.5 (dotted line). 
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